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Abstract 

Lysine acetylation is a frequently occurring posttranslational modi- 
fication; however, little is known about the origin and regulation 
of most sites. Here we used quantitative mass spectrometry to 
analyze acetylation dynamics and stoichiometry in Saccharomyces 
cereuisiae. We found that acetylation accumulated in growth- 
arrested cells in a manner that depended on acetyl-CoA generation 
in distinct subcellular compartments. Mitochondrial acetylation 
levels correlated with acetyl-CoA concentration in uiuo and 
acetyl-CoA acetylated lysine residues nonenzymatically in uitro. 
We developed a method to estimate acetylation stoichiometry and 
found that the vast majority of mitochondrial and cytoplasmic 
acetylation had a very low stoichiometry. However, mitochondrial 
acetylation occurred at a significantly higher basal level than 
cytoplasmic acetylation, consistent with the distinct acetylation 
dynamics and higher acetyl-CoA concentration in mitochondria. 
High stoichiometry acetylation occurred mostly on histones, 
proteins present in histone acetyltransferase and deacetylase 
complexes, and on transcription factors. These data show that 
a majority of acetylation occurs at very low levels in exponentially 
growing yeast and is uniformly affected by exposure to acetyl-CoA. 
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Introduction 

Lysine acetylation is an evolutionarily conserved, reversible post- 
translational modification (PTM] that is known to regulate protein 
function by site-specific acetylation (catalyzed by acetyltransferases) 
and deacetylation (catalyzed by deacetylases] . Acetyltransferases 
use acetyl-coenzyme A (acetyl-CoA) as a cofactor donating an acetyl 



group to the s-amino group of lysine. Acetylation is a well-known 
regulatory PTM in the context of nuclear signahng, in particular for 
regulating gene expression via modification of histones. The role of 
acetylation in regulating nuclear processes is consistent with the 
nuclear localization of most acetyltransferases and proteins with 
acetyllysine-binding bromodomains. A prominent role for non- 
nuclear acetylation was suggested by numerous proteomic studies 
showing that lysine acetylation occurs at thousands of sites through- 
out eukaryotic cells (Kim et al, 2006; Choudhary et al, 2009; Zhao 
et al, 2010; Weinert et al, 2011; Chen et al, 2012; Henriksen et al, 
2012; Lundby et al, 2012; Hebert et al, 2013). 

Notably, acetylation has been shown to occur frequently on 
mitochondrial proteins and a large number of sites are regulated 
by the sirtuin-class deacetylase-3 (SIRT3; Sol et al, 2012; Hebert 
et al, 2013; Rardin et al, 2013). These observations, and a number 
of supporting studies (Newman et al, 2012), have established acet- 
ylation as an important regulator of mitochondrial metabolism. 
However, the mechanisms by which acetylation occurs in this 
organelle are mostly unknown. Several recent reviews speculated 
that nonenzymatic acetylation by acetyl-CoA may be widespread 
in mitochondria (Guan & Xiong, 2010; Newman et al, 2012). Acet- 
ylation could regulate metabolism through a nonenzymatic mecha- 
nism, in which case determining the stoichiometry of acetylation is 
of critical importance for gauging the effects of acetylation on 
protein function. 

In this study we used quantitative mass spectrometry (MS) to 
measure changes in lysine acetylation abundance in Saccharomyces 
cereuisiae (budding yeast). We found that growth-arrest, combined 
with ongoing metabolism, resulted in the accumulation of acetyla- 
tion in a manner that depended on acetyl-CoA generation in distinct 
subcellular compartments. Acetylation dynamics in mitochondria 
correlated with acetyl-CoA levels in this compartment and we found 
that acetyl-CoA nonenzymatically acetylated protein in vitro. We 
developed a novel method to estimate acetylation stoichiometry and 
discovered that the vast majority of acetylation occurred at a very 
low level. These data provide the first global analysis of acetylation 
stoichiometry and indicate that metabolism regulates acetylation 
levels in distinct subcellular compartments through acetyl-CoA 
generation. 
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Results 

Experimental strategy for quantitative analysis of proteome and 
PTM dynamics 

We used stable isotope labeling with amino acids in cell culture 
(SILAC; Ong et al, 2002) to quantify differences in protein, acetyla- 
tion, and phosphorylation abundance by MS. Proteins from whole 
cell lysates were digested to peptides and acetylated peptides 
enriched using a polyclonal anti-acetyllysine antibody as previously 
described (Kim et al, 2006; Weinert et al, 2013a). Peptide fractions 
were analyzed by reversed-phase liquid chromatography coupled to 
high resolution liquid chromatography-tandem mass spectrometry 
[LC-MS/MS) and raw MS data were computationally processed 
using MaxQuant (Cox et al, 2011). The quantitative MS experiments 
performed in this study are summarized in Supplementary Table SI, 
and the raw data files have been deposited to the ProteomeXchange 
(Vizcaino et al, 2013), with the identifier PXD000507. 

Acetylation dynamics in growth-arrested yeast cells 

In previous work we showed that acetylation accumulates in 
growth-arrested bacteria (£. colt] cells due to prolonged exposure 
to, and a higher concentration of, acetyl-phosphate (Weinert et al, 
2013a). In order to test whether growth-arrest affects PTMs in a 
eukaryotic system we compared acetylation, phosphorylation, and 
protein levels between exponentially growing and stationary phase 
S. cerevisiae cells. We quantified more than 2500 acetylation sites 
(Supplementary Table S2), 3300 proteins (Supplementary 
Table S3), and 5000 phosphorylation sites (Supplementary 
Table S4) with a high correlation between biological replicates 
(Supplementary Figs SIA-C). Strikingly, stationary phase cells 
showed increased ( > 2-fold) acetylation at a majority (-70%) of 
quantified acetylation sites (median threefold increased; Fig 1). In 
contrast, protein and phosphorylation abundance, while affected in 
stationary phase cells, was not globally increased (Fig 1), indicat- 
ing that stationary phase did not cause the accumulation of 
proteins or PTMs generally. Furthermore, Gene Ontology (GO) 
enrichment analysis of protein and phosphorylation site changes 
revealed both up-regulation and down-regulation of specific pro- 
cesses in stationary phase cells (Supplementary Figs SID and E), 
suggesting that such changes occurred in a regulated manner. 
Using subcellular localization data from GFP-tagged proteins (Huh 
et al, 2003), we found that acetylation increased to the greatest 
degree on mitochondrial proteins (median 5.7-fold), which was 
significantly higher than the median increases seen for cytoplasmic 
(2.5-fold) and nuclear (1.5-fold) proteins (Fig 1). Increased acetyla- 
tion of mitochondrial proteins was nearly comprehensive (93% of 
sites were > 2-fold elevated; Fig 1 and Supplementary Table S2). 

During glycolysis, glucose is converted to pyruvate, which enters 
mitochondria and is further converted to acetyl-CoA by the pyruvate 
dehydrogenase complex (PDC), pyruvate can also be converted to 
acetyl-CoA in the cytoplasm via an acetate intermediate (Fig 2A) . In 
order to control the transition to stationary phase we growth- 
arrested exponentially growing cultures by transferring the cells into 
media containing different carbon sources, and lacking an essential 
amino acid (lysine), for -24 h. Notably, mitochondrial acetylation 



o 



o 



.2 Log2 SILAC ratio SP to EP 

? -4 -3 -2 -1 0 1 2 3 4 
E 



Proteins 3325 -0.1 1.0 

Acetyl proteins 959 -0.1 0.9 

Acetyl sites 2640 1.6 3.0 

Phosptio sites 6252 -0.1 0.9 t 

Mito. 498 2.5 5.7 

Cyto. 871 1.4 2,6 

Nuc. 265 0.5 1.5 



Acetyl sites 



•-HD-^ 
i-HD-i-i 

•-HH-H 
HZIZI-^^ 



p < 2e " 



p < 2e " 



I I I I I I I I — I 

^-3-2-101234 



Figure 1. Acetylation is globally increased in stationary phase yeast cells. 

Box plots showing the distributions of SII_AC ratios comparing stationary phase (SP) 
to exponential phase (EP) yeast for all quantified proteins (proteins), acetylated 
proteins (acetyl proteins), acetylation sites corrected for protein abundance 
changes (acetyl sites), and phosphorylation sites (phospho sites). The box portion of 
the plot indicates the middle 50% of the distribution, inner hatch marks denote 
9-91%, and whisker ends 2-98%, outliers are not shown. Acetylation sites 
occurring on proteins localized exclusively to mitochondria (Mito.), the cytoplasm 
(Cyto.), or the nucleus (Nuc.) are shown. Statistical significance was calculated 
using a Wilcoxon test, data is from two biological replicates. 



was unaffected when cells were growth-arrested in the presence of 
2-deoxy-D-glucose (2DG; a glucose analog that cannot be metabo- 
lized beyond the first step of glycolysis), indicating that increased 
mitochondrial acetylation depended on glycolysis (Fig 2B, Supple- 
mentary Fig S2A, and Supplementary Table S5). However, nuclear 
acetylation was globally reduced in the presence of 2DG (median 
3.5-fold), suggesting that nuclear deacetylases remain active under 
these conditions. 

In order to examine the dependence of mitochondrial protein 
acetylation on glycolysis, we quantified acetylation dynamics in 
pdalA cells that are unable to convert pyruvate to acetyl-CoA 
(Wenzel et al, 1992), or in citlA cells that are unable to convert 
acetyl-CoA to citrate (Kispal et al, 1988; Fig 2A). Mitochondrial 
acetylation was drastically reduced in exponentially growing pdalA 
cells (median SILAC ratio pda J zl/wild-type = 0.14) and was slightly 
increased in exponentially growing citlA cells (Fig 2C, Supplemen- 
tary Fig S2B, and Supplementary Table S5). Loss of Pdal com- 
pletely suppressed the increased acetylation of mitochondrial 
proteins in growth-arrested cells while loss of Citl further increased 
the acetylation of mitochondrial proteins (an additional 1.8-fold) in 
growth-arrested cells (Fig 2D, Supplementary Fig S2C, and Supple- 
mentary Table S7). Loss of Citl blocks the entry of acetyl-CoA into 
the citric acid cycle, thus, further increased acetylation in growth- 
arrested citl A cells is likely due to greater accumulation of acetyl- 
CoA in the mitochondria of these cells. We quantified half as many 
acetylation sites on mitochondrial proteins in pdal A cells compared 
to citl A cells, while the frequency of quantified sites on cytoplasmic 



2 Molecular Systems Biology 10 716 | 2014 



O 2014 The Authors 



Brian T. Weinert et al Acetylation dynamics and stoichiometry 



Molecular Systems Biology 



cytoplasm ,■' mitochondria 



acetyl-CoA , _ 
\ acetate - 

glycolysis | 

glucose pyruvate - 



' pyruvate acetyl-CoA — * citrate 

Pdal cm ■ ^ 



tWGA glucose 
wt GA 2DG 

wt GA glucose 
wt GA 2DG 

wt GA glucose 
wt GA 2DG 



216 
470 
142 
216 
470 
142 



c- 



54 
95 
388 



pda1A EP 

ciHA EP 
prfafA EP 

citIA EP 277 
pda1& EP 221 

citIA EP 151 



ivIGA 
pda1A GA 
Cit1A GA 
ivtGA 
pdalA GA 
citIA GA 
vWGA 
pda1A GA 
cit1A GA 



lucose 
lucose 
glucose 
glucose 
lucose 
lucose 
lucose 
lucose 
lucose 



E- 



wt GA glucose 
wt GA acetate 
wt GA glucose 
wt GA acetate 
wt GA glucose 
wt GA acetate 



308 
149 
274 
625 
544 
498 
300 
273 
231 

265 
223 
487 
660 
239 
316 



2 
o 
3 

w 



E 

2.5 
-0,1 
1.1 
-0.4 

0.0 
-1.8 

-2.8 
0.3 
-0.2 
0.0 
-0.1 
0.0 

2.6 
0.0 
3.3 
0.7 
1.0 
1.1 
0.2 
0.4 
0.4 

2.9 
1.1 
1.1 
1.7 
0.3 
1.0 



E 

5.6 
1.0 
2.2 
0.8 

1.0 
0.3 

0.1 
1.2 
0.9 
1.0 
0.9 
1.0 

6.0 
1.0 
9.6 
1.6 
2.1 
2.2 
1.1 
1.3 
1.3 

7.6 
2.2 
2.1 
3.2 
1.2 
2.0 



B Mitocliondhal 
I I Cytoplasmic 
I Nuclear 

Median Log2 SILAC ratio 
-3-2-10123 



Mito. 
Milo. 
Cyto. 
Cyto. 
Nuc 
Nuc 

Mito 
Mito 
Cyto. 
Cyto. 
Nuc, 
Nuc, 

Mito, 
Mito, 
Mito, 
Cyto, 
Cyto, 
Cyto, 
Nuc, 
Nuc 
Nuc 

Mito, 
Mito, 
Cyto, 
Cyto, 
Nuc, 
Nuc 




Figure 2. Quantitative analysis of acetylation dynamics in yeast. 

A Model showing the formation of acetyl-CoA from glucose and acetate, 
key enzymes are shown in red type, 

B-E The figure shows the conditions analyzed in each experiment, the cell 
type (wild-type (lut) or indicated mutant strains), the growth state 
[exponential phase (EP) or growth-arrested (GA)], the number of 
acetylation sites analyzed (# sites), the median Log2 and linear SILAC 
ratios comparing the indicated condition to wild-type EP cells, and the 
subcellular localization of the analyzed acetylation sites on proteins 
localized to mitochondria (Mito.), the cytoplasm (Cyto.), or the nucleus 
(Nuc). Cells were growth-arrested by transferring an exponential phase 
culture into media lacking lysine and containing the indicated carbon 
sources, glucose, acetate, or 2-deoxy-D-glucose (2DG). The bar chart 
shows the median Log2 SILAC ratios comparing the indicated condition 
to wild-type EP cells, statistical significance was calculated by Wilcoxon 
test. Increased acetylation requires glycolysis (B). Data is from two 
biological replicates. Mitochondrial acetylation in exponentially growing 
cells requires Pdal (C). Increased mitochondrial acetylation in 
growth-arrested cells is suppressed by loss of Pdal and enhanced by 
loss of Citl (D). Data from two biological replicates is shown, significant 
differences are relative to wild-type cells. Acetate promotes cytoplasmic 
and nuclear acetylation (E). Data is from two biological replicates. 



and nuclear proteins was not altered (Figs 2C and D], indicating 
that acetylation of many sites in pdal A cells had diminished to 
undetectable levels in mitochondria specifically. 

We next compared growth-arrested yeast in the presence of glu- 
cose or acetate to test whether acetate would alter the acetylation 
dynamics in growth-arrested cells. Consistent with conversion of 
acetate to acetyl-CoA by acetyl-CoA synthetase 2 (Acs2; Fig 2A; 
Takahashi et al, 2006) in the cytoplasm, acetate preferentially 
promoted the acetylation of cytoplasmic and nuclear proteins 
(Fig 2E, Supplementary Fig S2D, and Supplementary Table S8). 

Acetylation levels in mitochondria correlate with 
acetyl-CoA concentration 

We showed that mitochondrial acetylation was suppressed by dele- 
tion oipdal, suggesting that most mitochondrial acetyl-CoA is gener- 
ated in a Pdal -dependent manner in the growth conditions used in 
our study. We compared acetyl-CoA levels between wild-type, 
pdalA, and citlA cells (Fig 3A}. Acetyl-CoA was reduced by approxi- 
mately 1/3 in pdal A cells, indicating that the mitochondrial 
acetyl-CoA pool constitutes at least 1/3 of the total acetyl-CoA pool. 
However, mitochondria occupy -1-2% of the total cellular volume in 
yeast (Uchida et al, 2011), indicating that acetyl-CoA levels are 
substantially (~20-30-fold) higher in mitochondria compared to the 
cytoplasm and nucleus. We hypothesized that increased acetylation 
in growth-arrested cells may occur due to higher concentrations of 
acetyl-CoA or due to prolonged exposure of proteins to acetyl-CoA. 
We found that acetyl-CoA levels were reduced in growth-arrested 
cells (Fig 3A), suggesting that acetylation occurs due to prolonged 
exposure to acetyl-CoA. However, this difference may be due to 
reduced recovery of acetyl-CoA from growth-arrested cells, which 
have a substantially increased cell wall and are known to be 
refractory to cell lysis. We consistently recovered less protein from 
growth-arrested cells (data not shown), suggesting that these cells 
were more difficult to lyse, or that cell size and/or protein content 
was reduced under these conditions. Such differences in cell physiol- 
ogy may explain the lower amount of acetyl-CoA per ODgoo of 
growth-arrested cells as determined by our assay. Regardless, acetyl- 
CoA levels were elevated in growth-arrested citl A cells, consistent 
with our observation of increased acetylation in these cells. If we 
assume that the reduction of acetyl-CoA in pdalA cells indicates the 
mitochondrial acetyl-CoA pool in these cells, then mitochondrial 
acetyl-CoA in citl A cells was increased ~1. 9-fold and ~1. 7-fold after 8 
and 16 h growth-arrest, respectively (Fig 3A), mirroring the 1.8-fold 
increase in mitochondrial acetylation in growth-arrested citlA cells 
(Fig 2D). 

The physiological concentration of acetyl-CoA in yeast is not 
well-studied. One study estimated acetyl-CoA levels (3-30 |.tM) in 
nutrient-starved yeast undergoing metaboUc cycles (Cai et al, 2011), 
conditions that contrast with the excess of glucose and high growth 
rates of yeast grown on synthetic complete (SC) media in our study. 
We estimated a similar cellular concentration of -30 |j.M acetyl-CoA 
in exponentially growing cells. However, this estimate assumes 
complete recovery of acetyl-CoA, and is therefore likely to under- 
estimate the actual cellular concentration. We showed that 
mitochondrial acetyl-CoA was -20-30-fold higher than non- 
mitochondrial acetyl-CoA, suggesting a concentration of acetyl-CoA 
in mitochondria that approaches the millimolar range (-0.5-1 mM 
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Figure 3. Acetyl-CoA concentration in cells and nonenzymatic acetylation by acetyl-CoA. 

A Acetyl-CoA concentration was determined in the indicated cell types during exponential phase (EP) growth or after the indicated time of growth-arrest (GA) in the 

presence of glucose. Data are from two independent biological replicates. 
B The bar graph shows the abundance of the indicated, acetylated (ac) peptides relative to an untreated control sample. Error bars indicate standard deviation of the 

indicated number (n) of independently quantified peptides. The significance (P) of increased acetylation at 100 |.iM acetyl-CoA was calculated by two-tailed t-test 

assuming equal variance. 

C The column graph shows the relative abundance of non-acetylated peptides. Error bars indicate standard deviation of the indicated number (n) of independently 
quantified peptides. 



based on our estimates]. This estimate is consistent with previous 
work showing that acetyl-CoA can reach millimolar levels in rat 
liver mitochondria (Garland et al, 1965). 

Acetyl-CoA nonenzymatically acetylates protein 

A study in 1970 showed that acetyl-CoA could nonenzymatically 
acetylate lysine rich histone fractions and synthetic polylysine at a 
concentration of -0.5 mM, with increased kinetics at higher pH, and 
an activation energy of 7.5 kcal, suggesting that the reaction could 
occur under physiological conditions (Paik et al, 1970). We sought 
to independently verify nonenzymatic acetylation by acetyl-CoA 
using mass spectrometry and a quantitative strategy based on iso- 
baric mass tags (Supplementary Fig S3). Bovine serum albumin 
(BSA) was treated with increasing concentrations of acetyl-CoA at 
physiological pH, modestly, but significantly, increased acetylation 
was seen on several peptides at a concentration of 100 |.tM acetyl- 
CoA, while 1 and 10 mM acetyl-CoA caused increased acetylation 
on all peptides analyzed (Fig 38). In contrast, the abundance of 
non-acetylated peptides was unaffected by acetyl-CoA (Fig 3C) . This 
acetylation is likely to be nonenzymatic since we treated a purified 
extracellular serum protein (BSA) and acetylation did not occur 
until acetyl-CoA concentration was well in excess of the known 
range of binding affinities of acetyltransferases for acetyl-CoA (Alb- 
augh et al, 2011a). 

Acetylation stoichiometry in yeast 

The striking increase (median 9.5-fold) in mitochondrial acetylation 
in growth-arrested citlA cells indicated that mitochondrial 
acetylation occurred at a low-level (maximum median stoichiometry 
of -10% in exponentially growing cells). One method used to 
analyze PTM stoichiometry on a global scale is to compare the 



relative abundances of posttranslationally modified and unmodified 
corresponding peptides (CPs; Olsen et al, 2010; Wu et al, 2011). If 
the abundance of a modified peptide is substantially altered then the 
abundance of the CP should be affected in an inverse manner. How- 
ever, CP abundance will only be measurably affected if PTM stoichi- 
ometry is relatively high. Since our results suggested a low 
stoichiometry of modification, we devised a strategy to assay acety- 
lation stoichiometry based on partial chemical acetylation of lysine 
residues. Partial chemical acetylation will cause the greatest relative 
increases in acetylation on sites with the lowest initial stoichiome- 
try, while sites with higher initial stoichiometry will be partially or 
not at all affected. In order to find conditions that resulted in partial 
chemical acetylation, we chemically acetylated bovine serum albu- 
min (BSA) with acetyl-phosphate (AcP) or acetic anhydride. We 
found that AcP caused partial chemical acetylation that did not mea- 
surably affect the abundance of CPs, while acetic anhydride caused 
comprehensive acetylation and resulted in dramatically reduced 
abundance of CPs (Supplementary Fig S4 and supplemental text) . 

In order to examine acetylation site stoichiometry in S. cerevisiae, 
we used a SILAC-based quantitative approach to quantify chemi- 
cally acetylated and untreated peptides in the same MS experiment 
(Fig 4). An untreated control sample was prepared in order to iden- 
tify acetylation sites that occurred in the absence of AcP-treatment 
and to determine the quantitative variability in our experiments. 
Treatment with AcP resulted in substantially increased (>2-fold) 
acetylation at 74% (10 mM AcP) and 88% (100 mM AcP) of acety- 
lation sites (Fig SA and Supplementary Table S9). In contrast, <1% 
of acetylation sites were substantially increased in the untreated 
(0 mM AcP) lysate, reflecting the quantitative accuracy of our 
assays (Fig 5A). An independent experimental replicate was per- 
formed with untreated (0 mM) and 100 mM AcP-treated lysate 
(Supplementary Table S9). Acetylation changes were similar in 
experimental repUcates, as shown by a Spearman's correlation 
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Figure 4. Method used to assay acetyl-phosphate sensitivity in 
yeast lysate. 

Yeast proteins were metabolically labeled with unlabeled "light" lysine or with 
"heavy" isotope labeled lysine. Protein lysates were treated with acetyl- 
phosphate (AcP) or were mock-treated by addition of H2O. Equal amounts of 
protein were mixed and digested to peptides with trypsin protease. Acetylated 
(Ac) peptides were immune-affinity enriched using agarose-coupled anti- 
acetyllysine polyclonal antibody and analyzed by MS. Increased acetylation 
from AcP treatment causes an increase in the relative abundance of the light 
labeled peptide, enabling quantification of acetylation changes induced by AcP. 



coefficient of 0.89 [Supplementary Fig S5). We identified 9415 acet- 
ylation sites in 100 mM AcP-treated lysates, 1765 sites had SILAC 
ratios, and 916 of tiiese sites were previously identified (Supplemen- 
tary Table S2), or were found in untreated control samples, indicat- 
ing that these sites are naturally occurring acetylation sites in yeast. 

Using only the 916 naturally occurring sites with SILAC ratios, 
we found that 65% of the sites were more than 10-fold increased in 
acetylation (Fig 5B), indicating maximum initial stoichiometry that 
was < 10% at these sites. However, acetylation stoichiometry at 
these positions is hkely to be lower than 10% (which assumes com- 
prehensive acetylation) because our analysis of AcP-treated BSA 
(Supplementary Fig S4], showed that 100 mM AcP caused partial 
chemical acetylation without a reduction in the abundance of CPs. 
CP abundance was similarly unaffected in 100 mM AcP-treated 
yeast lysate (Supplementary Fig S6], indicating that AcP-treatment 
caused a low-level of partial chemical acetylation. 

In order to better characterize the degree of partial chemical acet- 
ylation by 100 mM AcP, we used an absolute quantification (AQUA) 
method (Gerber et al, 2003), to quantify acetylation levels. Acetylated 
and unmodified peptide intensities were compared to heavy-labeled 
peptide standards, allowing us to determine acetylation stoichiometry 
directly. We analyzed eight acetylation sites with varying sensitivity 
to 100 mM AcP on two proteins, 3-phosphoglycerate kinase (Pgkl) 



and fatty acid synthetase (Fas2; Table 1). Surprisingly, we were 
unable to detect acetylated peptides from purified Pgkl, and found 
just one acetylated peptide from purified Fas2, even though we 
detected unmodified peptides covering 89 and 72 % of the Pgkl and 
Fas2 protein sequences, respectively. In order to detect acetylation 
on these proteins we treated cell lysates with 100 mM AcP and puri- 
fied Pgkl and Fas2 for AQUA analysis. Comparison to heavy-labeled 
peptide standards indicated that acetylation stoichiometry, after 
AcP-treatment, was < 1 % at all eight sites, with a median degree of 
chemical acetylation that was just 0.07% (Table 1, Supplementary 
Figs S7A-D). AcP-sensitivity was well-correlated with acetylation 
stoichiometry as determined by the AQUA method (Spearman's cor- 
relation of —0.92, Fig 5C), confirming our prediction that low stoi- 
chiometry sites would be most sensitive to partial chemical 
acetylation and providing independent validation of our method. 

We estimated acetylation stoichiometry based on the conserva- 
tive assumption that chemical acetylation from 100 mM AcP was 

< 1% at all sites. A site with 10-fold increased acetylation after 
AcP-treatment was estimated to have a stoichiometry that is < 0.1% 
while a site with 20-fold increased acetylation was estimated to have 
a stoichiometry that is < 0.05%. Using this approach we estimated 
acetylation stoichiometry based on the observed SILAC L/H ratios 
after treatment with 100 mM AcP and found that 25% of sites were 

< 0.02% acetylated, 50% of sites were < 0.05% acetylated, 66% of 
sites were < 0.1% acetylated, and a remarkable 86% of sites were 

< 1 % acetylated (Fig 5D) . These estimates contrast with previously 
determined phosphorylation site stoichiometries in yeast, where 
89% of sites were estimated to have stoichiometries that were >1 % 
(Wu et al, 2011; Fig 5E), indicating that a much greater fraction of 
acetylation occurs with a low stoichiometry. It was not possible to 
accurately estimate the stoichiometry of sites that were AcP-insensi- 
tive (SILAC ratio L/H < 2) as the relative changes in acetylation 
were of a similar magnitude to the variability of these measure- 
ments. Thus, AcP-insensitive sites were estimated to have acetyla- 
tion stoichiometry that is >1 % . 

In order identify independent parameters that could confirm our 
stoichiometry estimates we compared acetylated peptide intensity 
and abundance-corrected acetylated peptide intensity to AcP-sensi- 
tivity (Supplementary Fig S8). We found that correcting acetylated 
peptide intensity with protein abundances determined by an inten- 
sity-based absolute quantification (IBAQ) method (Schwanhausser 
et al, 2011), could best distinguish AcP-insensitive sites (SILAC ratio 
L/H < 2) from AcP-sensitive sites (SILAC ratio L/H > 2; Fig 5F and 
Supplementary Fig S8) and had the best overall correlation 
(p = —0.57) with AcP sensitivity (Supplementary Fig S8). These 
data independently confirmed our stoichiometry estimates based on 
AcP-sensitivity, as this parameter was calculated using peptide 
intensities and protein abundances determined in untreated control 
samples. 

We estimated acetylation stoichiometry for the 916 sites with 
SILAC ratios; however, an additional 1540 sites without SILAC 
ratios were independently identified as naturally occurring sites 
(identified in untreated control samples or in Supplementary 
Table S2). SILAC ratios were estimated for these sites by calculat- 
ing the increased intensity of AcP-treated "light" peptides relative 
to an empirically determined detection limit for naturally occurring 
"heavy" peptides (see Materials and Methods). Using this method 
we calculated minimum ratios of increased acetylation at 1516 
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Figure 5. Most acetylation sites are modified withi a low stoichiometry. 

A The majority of yeast acetylation sites are highly sensitive to partial chemical acetylation by AcP. The histogram shows the distribution of SILAC lyH ratios for the 
indicated samples. 

B AcP caused substantially increased acetylation at a majority of sites. The histogram shows acetylation changes induced by 100 mM AcP in two experimental 

replicates, only sites that were independently identified in cells without AcP treatment are shown. 
C Acetylation stoichiometry is inversely proportional to AcP-sensitivity. The scatterplot shows the relationship between AcP-sensitivity (SILAC ratio L/H 100 mM AcP) 

and acetylation stoichiometry (LoglO stoichiometry) as determined by AQUA analysis (Table 1). The Spearmans correlation (p) and the significance by two-tailed 

test (P-value) are shown. 

D Most acetylation occurs with a low stoichiometry. Absolute acetylation site stoichiometries were estimated based on relative abundance changes (SILAC ratio L/H) 

after treatment with 100 mM AcP and using an estimate that AcP treatment caused < 1% chemical acetylation. 
E For comparison, previously determined phosphorylation site stoichiometries are shown (Wu et al, 2011). 

F iBAQ-based abundance corrected acetylated peptide intensity (l/iBAQ-A) is proportional to AcP sensitivity. The box plots show the distributions of l/iBAQ-A values for 

the indicated classes of acetylation sites. AcP-insensitive sites had a significantly (p) higher distribution of l/iBAQ-A values compared to either AcP-sensitive sites or 

sites without SILAC ratios. Significance was calculated by Wilcoxen test 
C Sites without SILAC ratios are highly sensitive to AcP. The minimum ratio of increased acetylation was determined by calculating the increased intensity of AcP- 

treated "light" peptides relative to an empirically determined detection limit for "heavy" SILAC peptides (see Materials and Methods). 
H Absolute acetylation stoichiometry of sites without SILAC ratios was estimated to be very low. Stoichiometry was estimated by the same method used to estimate 

stoichiometry of sites with SILAC ratios in (D) using the minimum ratios of increased acetylation shown in (C). 



sites and found that AcP caused a median 70-fold increase in acet- 
ylation (Fig 5G). Assuming that 100 mM AcP caused < 1% chemi- 
cal acetylation, the substantially increased acetylation at these 
sites indicated a very low stoichiometry of acetylation (Fig 5H) 
and this finding was independently verified by examining the 



abundance-corrected peptide intensity (I/iBAQ-A) distributions for 
these sites (Fig 5F). By including these data we estimated acetyla- 
tion stoichiometry at a total of 2432 sites in exponentially growing 
yeast and found that 95% of all acetylation occurs at an estimated 
stoichiometry that is < 1 % . 
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Table 1. Absolute quantification (AQUA) of acetylated peptides from Pgl<l and Fas2. Acetylation stoichiometry was determined by AQUA 
for Pgkl and Fas2 after treatment with 100 mM AcP (see Supplementary Fig S6). Initial stoichiometry was calculated by dividing the 
stoichiometry in 100 mM AcP by the Ratio L/H 100 mM AcP and the degree of chemical acetylation is the difference between the initial 
stoichiometry and the stoichiometry in 100 mM AcP 



Protein 


Peptide sequence 


Stoichiometry 

inn tvtKA A^D lo/\ 

lUO mM ACr \ /o) 


Ratio L/H 
lUU mjvl ACr 


Initial 

stoichiometry (%) 


Degree chemical 
acetylation (%) 


Pgkl 


AACFLLEK(ac)ELK 


0.023 


80.8 


0.0003 


0.02 


Pgkl 


ACAEIVPK(ac)LMEK 


0.095 


41.5 


0.0023 


0.09 


Pgkl 


FAACTK(ac)ALLDEWK 


0.119 


27.7 


0.0043 


0.11 


Fas2 


QVLDVDPVYKDVA(ac)PTCPK 


0.032 


108.3 


0.0003 


0.03 


Fas2 


LIEPELFNCYNPE(ac)K 


0.092 


33.5 


0.0027 


0.09 


Fas2 


GATLYIP{ac)ALR 


0.068 


11.3 


0.0060 


0.06 


Fas2 


SECNPVIGVFQ(ac)FLTCHPK 


0.016 


5.3 


0.0031 


0.01 


Fas2 


AND(ac)NESATINEMMK 


0.181 


1.7 


0.1055 


0.08 



Functional analysis of proteins with high stoichiometry 
acetylation sites 

In order to characterize proteins with AcP-insensitive, high stoichi- 
ometry sites we used GO term enrichment analysis. Proteins with 
AcP-insensitive sites were significantly more frequently associated 
with nuclear processes (Fig 6A). More than half of all AcP-insensi- 
tive sites occurred on histones, proteins present in histone deacety- 
lase (HDAC) or histone acetyltransferase (HAT] complexes, and on 
transcription factors. We compared the effect of AcP-treatment on 
proteins present in different subcellular compartments [Fig 6B). 
Consistent with GO term analysis, 77% of the high-stoichiometry, 
AcP-insensitive sites occurred on nuclear proteins. In contrast, 97% 
of mitochondrial and 94% of cytoplasmic acetylation sites were 
AcP-sensitive, indicating that the vast majority of acetylation sites in 
these subcellular compartments were acetylated at a low stoichiom- 
etry. Notably, mitochondrial acetylation sites were significantly 
(p = 5.8e"", Wilcoxon test) less affected (median 14-fold increased] 
by AcP than cytoplasmic sites [median 30-fold increased; Fig 6B}, 
suggesting a higher basal level of acetylation in mitochondria. This 
finding was independently confirmed by a significantly higher distri- 
bution of acetylated peptide I/iBAQ-A for mitochondrial proteins 
than cytoplasmic proteins (Fig 6C). In contrast, non-acetylated pep- 
tide I/iBAQ-A was similarly distributed between mitochondrial and 
cytoplasmic proteins (Fig 6C) . It is also notable that the distribution 
of AcP-sensitivity for nuclear proteins was bi-modal, indicating 
distinct populations of high and low stoichiometry acetylation sites 
in the nucleus. The subcellular localization of high-stoichiometry, 
AcP-insensitive sites was significantly different to low-stoichiome- 
try, AcP-sensitive sites and low-stoichiometry sites without SILAC 
ratios in all three compartments (Fig 6D) . 

MS analysis is biased towards abundant proteins, which are 
more readily detected in the MS. Nearly every abundant protein 
identified by MS was also acetylated [Fig 6E), and the frequency of 
detected acetylation sites was proportional to protein abundance 
[Supplementary Fig S9A], indicating that protein abundance is a 
limiting factor in the identification of many sites. The technical bias 
to detect acetylation on abundant proteins was less pronounced for 
mitochondrial proteins compared with cytoplasmic proteins, consis- 
tent with our estimation of higher acetylation levels in mitochondria 



(Fig 6E}. Similarly, nuclear proteins were least biased to detect acet- 
ylation on abundant proteins, as this subcellular compartment con- 
tained the most high stoichiometry acetylation sites (Fig 6E). Thus, 
the abundance of proteins with observed acetylation sites is consis- 
tent with our stoichiometry estimates. 

Acetylated sites with high stoichiometry are likely to be impor- 
tant for protein function. Many previously known regulatory acety- 
lation sites were AcP-insensitive, including Smc3 lysines 112 and 
113 (Zhang et al, 2008), Sas2 (MYST homolog) lysine 168 (Yuan 
et al, 2012], Yng2 lysine 170 (Lin et al, 2008), RTT109 lysine 290 
(Albaugh etal, 2011b), SNF2 lysines 1494 and 1498 [Kim et al, 
2010), histone H2AZ (Htzl) lysines 4, 9, 11, and 15 [Babiarz et al, 
2006; Millar et al, 2006; Lin et al, 2008), histone H3 (Hhtl) lysines 
19, 24, and 57 (Zhang et al, 1998; Suka et al, 2001; Hyland et al, 

2005) , histone H4 (Hhfl) lysines 6, 9, 13, and 17 (Suka et al, 2001) 
and histone H2B [Htb2) lysines 16 and 17. These 20 sites constitute 
18% of the 111 AcP-insensitive sites that we identified, indicating 
that AcP-insensitivity is a good predictor of functionally important 
acetylation sites. The enrichment of functionally characterized sites 
within the group of AcP-insensitive sites is highly significant 
(p = le"'^. Fisher exact test). Furthermore, known functional sites 
had significantly higher I/iBAQ-A values, indicating that this param- 
eter is also useful in distinguishing functionally-relevant, high- 
stoichiometry sites (Supplementary Fig S9B and Supplementary 
Table S9). 

Discussion 

while acetylation has been identified on thousands of proteins in 
diverse subcellular compartments, the mechanisms of non-nuclear 
acetylation and its regulation are not well understood. We found 
that acetylation levels were dynamically affected by growth-arrest 
and the generation of acetyl-CoA in distinct subcellular compart- 
ments. Acetyl-CoA was previously shown to exist in distinct mito- 
chondrial and non-mitochondrial pools in yeast (Takahashi et al, 

2006) . We showed that mitochondrial acetylation occurs within 
mitochondria and is uniformly affected by growth-arrest and the 
generation of acetyl-CoA by the PDC. We found that acetyl-CoA 
levels correlated with acetylation changes in mitochondria and we 
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Figure 6. Functional analysis of liigli stoicliiometry acetylation sites. 

A AcP-insensitive acetylation sites occur on proteins associated with nuclear processes. Gene Ontology (CO) term enrichment was performed by comparing proteins 
with AcP-insensitive acetylation sites (ratio IVH < 2) to all acetylated proteins. The bar graph shows the percentage of AcP-insensitive sites occurring on proteins 
associated with the indicated GO terms. P-values (P) indicate the statistical significance of GO term enrichment by Fishers exact test 

B AcP-insensitive acetylation sites occur on nuclear proteins. The histogram shows the distribution of SILAC L/H ratios occurring on proteins localized to the indicated 
subcellular compartments. 

C Acetylated peptides from mitochondrial proteins have a significantly higher median l/iBAQ-A value compared to acetylated peptides from cytoplasmic proteins. The 
box plots show the l/iBAQ-A distributions for the indicated classes of peptides occurring on mitochondrial (Mito.) or cytoplasmic (Cyto.) proteins. Signficance (p) was 
determined by Wilcoxon test 

D AcP-insensitive sites (SILAC Ratio < 2) have a significantly different subcellular distribution. Sites without SILAC ratios (No ratio) have a similar subcellular 
distribution to AcP-sensitive sites (SILAC Ratio > 2). The bar graph shows the fraction of sites localized to the indicated subcellular compartments; mitochondria 
(Mito.), cytoplasm (Cyto.), or nucleus (Nuc). Significance (P) was calculated by Fisher exact test. 

E Detection of acetylation sites is biased to occur on abundant proteins and this bias is more pronounced for sites with the lowest estimated stoichiometries. The 
histograms show the distributions of iBAQ protein abundances for observed proteins (n = 3,104). The distributions of the indicated classes of acetylated proteins 
occurring exclusively in the indicated subcellular compartments is shown in red. The numbers of acetylated proteins are shown in parenthesis and the median 
LoglO iBAQ abundance for these acetylated proteins is shown. 



confirmed that acetyl-CoA can nonenzymatically acetylate proteins. 
Our analysis of acetylation stoichiometry showed that the vast 
majority of acetylation occurs at very low levels, an important 
observation with widespread implications for understanding non- 
nuclear acetylation. 

A few recent reviews speculated that acetylation may occur 
nonenzymatically by exposure to acetyl-CoA, particularly in mito- 
chondria where a protein acetyltransferase activity has not been 
definitively identified (Guan & Xiong, 2010; Newman et al, 2012]. 
Acetyl-CoA has been previously shown to nonenzymatically acety- 
late proteins (Paik et al, 1970] and the elevated pH in the mito- 
chondrial matrix would favor this reaction (Paik et al, 1970; Llopis 
et al, 1998; Abad et al, 2004; Wagner & Payne, 2013]. In addition, 
recent work showed that lysines can be nonenzymatically modified 
by other glycolytic intermediates, such as 1,3-bisphosphoglycerate 
(Moellering & Cravatt, 2013], and succinyl-CoA (Wagner & Payne, 
2013; Weinert et al, 2013b]. We found that nearly all acetylation 
occurred at a low-level and was uniformly regulated by manipula- 
tions that affected the generation of acetyl-CoA in distinct subcellu- 



lar compartments. Mutations that blocked acetyl-CoA generation or 
its downstream conversion to citrate had opposing effects on mito- 
chondrial acetylation in growth-arrested yeast cells, linking acetyl- 
CoA levels to mitochondrial acetylation levels. We furthermore 
showed that acetyl-CoA levels correlated with acetylation abun- 
dance in mitochondria. Growth-arrest in the presence of glucose or 
acetate affected acetylation levels in a manner that was consistent 
with separate mitochondrial and non-mitochondrial acetyl-CoA 
pools and generation. Mitochondrial proteins had a significantly 
higher basal level of acetylation than cytoplasmic proteins, consis- 
tent with our results indicating a higher concentration of acetyl- 
CoA in this organelle. Acetylation is widespread in yeast (Henrik- 
sen et al, 2012], yet the number of known acetyltransferases in 
yeast is small (-15; Cherry et al, 2012], particularly compared to 
the number of kinases catalyzing protein phosphorylation (-130; 
Manning et al, 2002; Cherry et al, 2012]. Moreover, most acety- 
ltransferases have known nuclear functions, consistent with our 
observation that the majority of high stoichiometry acetylation 
occurred in the nucleus. Thus, the global acetylation dynamics that 



8 Molecular Systems Biology 10 716 | 2014 



O 2014 The Authors 



Brian T. Weinert et al Acetylation dynamics and stoichiometry 



Molecular Systems Biology 



we observed in growth-arrested yeast, their dependence on distinct 
acetyl-CoA pools and generation, a preponderance of low stoichi- 
ometry acetylation, and the restriction of high stoichiometry acety- 
lation to the nucleus, are highly consistent with a nonenzymatic 
mechanism of acetylation. These observations contrast with the 
canonical model of protein regulation by site-specific PTMs and 
indicate that acetylation may act in a global manner to regulate 
proteins, or conversely, that most acetylation occurs as a low-level 
protein lesion that accumulates under specific conditions. How- 
ever, widespread low-level acetylation could occur due to a pro- 
miscuous acetyltransferase activity that is sensitive to changes in 
acetyl-CoA levels. 

To better understand our estimates of acetylation stoichiometry, 
it is important to note that individual lysines may vary in the degree 
of chemical acetylation by AcP. Without knowing the exact degree 
of chemical acetylation at each position it is not possible to make a 
precise estimate of acetylation stoichiometry. In order to account for 
this variability we based our estimates on the conservative assump- 
tion that the degree of partial chemical acetylation by AcP was <1 % . 
Since we determined that 100 mM AcP caused a median chemical 
acetylation of just 0.07% (between 0.01 and 0.11%), we are likely 
to overestimate the stoichiometry of acetylation. Thus, stoichiome- 
try estimates are presented as less than (<) values, to designate an 
estimated stoichiometry that is less than the indicated amount. In 
addition, some lysines may be inaccessible or otherwise unreactive, 
resulting in an incorrect estimate of high stoichiometry acetylation. 
Regardless of these limitations, the results obtained using our 
approach indicated that it was able to distinguish acetylation stoichi- 
ometry at the site level. Absolute quantification indicated that acety- 
lation stoichiometry was inversely proportional to AcP-sensitivity 
(Table 1, Fig 5C). The highly significant bias to identify high stoichi- 
ometry sites on nuclear proteins indicated that the identification of 
such sites was non-random. Similarly, sites we predicted to have 
high stoichiometry were significantly more likely to have known 
functional roles. Our method further distinguished between low 
stoichiometry acetylation in the cytoplasm and mitochondria, indi- 
cating that the difference in acetylation levels between these two 
subcellular compartments was greater than the inherent variability 
of lysine reactivity with AcP. Thus, the resolution of our assay was 
not limited by the variability in lysine reactivity. 

Phosphorylation is known to regulate proteins in a site-selective 
manner and is one of the most extensively studied and widespread 
PTMs in eukaryotic cells. The number of detected phosphorylation 
sites increases with organism complexity, reflecting its role in regu- 
lating diverse cellular processes. In contrast, acetylation occurs as 
frequently in bacteria as in human cells (Weinert et al, 2013a). Simi- 
lar to bacteria, acetylation is overrepresented in mitochondria com- 
pared with phosphorylation (Weinert et al, 2011). Thus the 
frequency of acetylation is uncoupled from organism complexity, 
suggesting that nonenzymatic acetylation may play a large role in 
the generation of these sites. We showed that basal acetylation lev- 
els are elevated in mitochondria, likely due to a higher concentra- 
tion of acetyl-CoA in this organelle. These observations suggest that 
the endosymbiotic evolution of mitochondria as the metabolic cen- 
ters of eukaryotic cells may have limited widespread acetylation 
outside of mitochondria by compartmentalizing acetyl-CoA gener- 
ated during metabolism, thereby enabhng the evolution of acetyla- 
tion signaling in the nucleus. 



Materials and Methods 

Yeast growth and lysate preparation 

Wild-type Saccharomyces cerevisiae cells (BY4742, MATalpha 
hisSAl leu2A0 lys2A0 ura3A0; ThermoFisher Scientific, Slangerup, 
Denmark), pdai /I (BY4742, pdaJ::KanMX; ThermoFisher Scientific), 
and citlA (BY4742, citJ::KanMX; Open Biosystems) were cultured in 
synthetic complete media (US Biological, Salem, MA, USA) supple- 
mented with '^Cs'^Ni-lysine (SILAC "light") or "Cs'^Nj-lysine 
(SILAC "heavy"). Cells were harvested at the indicated time points, 
washed once with sterile H2O, and resuspended in lysis buffer 
(50 mM Tris, pH7.5, 150 mM NaCl, 1 mM EDTA, Ix mini complete 
protease inhibitor cocktail (Roche, Basel, Switzerland), 5 mM 
sodium fluoride, 1 mM sodium orthovanadate, 5 mM beta-glycero- 
phosphate, 10 mM nicotinamide, and 5 ).iM tricostatin A) at 
-50 OD|5oo cells/ml lysis buffer. The cell suspension was frozen 
drop-wise in liquid nitrogen and ground in a liquid nitrogen 
chilled steel container by the Retsch MM 400 Ball Mill (Retsch, 
Haan, Germany) for 5 min at 25 Hz. The lysate was thawed, NP-40 
and sodium deoxycholate were added to a final concentration of 1 
and 0.1%, respectively, and clarified by centrifugation. The lysate 
supernatent was precipitated with four volumes — 20°C acetone. The 
acetone precipitate was dissolved in urea solution (6 M urea, 2 M 
thio-urea, 10 mM Hepes pH8.0) and protein concentration 
determined by Quick-Start Bradford assay (Bio-Rad, Copenhagen, 
Denmark) . 

Chemical acetylation of BSA and yeast lysate 

For treatment of BSA with acetyl-CoA, BSA was prepared in PBS 
(pH 7.2) at 10 mg/ml and mixed with 1/10 volume acetyl-CoA 
(Sigma, Copenhagen, Denmark) prepared fresh in H2O. Reactions 
were incubated at 30°C for 4 h and stopped by addition of nine vol- 
umes — 20°C acetone. Acetone precipitated protein was resuspended 
in 200 |il of 100 mM triethyl ammonium bicarbonate (TEAB) and 
digested to peptides by addition of 1/100 (w/w) trypsin protease 
(Sigma) and incubation at 37°C for 16 h. Tryptic peptides were 
labeled with TMT sixplex isobaric mass tags (ThermoFisher 
Scientific) according to the manufacturer's instructions. Chemical 
acetylation with acetic anhydride was performed as described (Guan 
et al, 2010). Briefly, 20 ml of 1 mg/ml BSA (Sigma) was prepared in 
0.1 M Na2C03. 100 |il acetic anhydride (Sigma) was added slowly 
over 10 min in a glass beaker with constant stirring, followed by 
400 ).d pyridine (Sigma) over 30 min, also with constant stirring. 
The reaction was then left to incubate at room temperature for 4 h 
and quenched by addition of 400 )il 1 M Tris-base (Sigma). Organic 
solvents were removed by five rounds of centrifuge filtration in an 
Amicon Ultra 15 (Millipore, Billerica, MA, USA) with a 10 kD 
molecular weight cutoff. For treatment with AcP, BSA was prepared 
in yeast lysis buffer (as above, with NP-40 and deoxycholate added) 
at a concentration of 20 mg/ml. Yeast lysate was similarly prepared 
at a concentration of 20 mg/ml. BSA or yeast lysate was then mixed 
with 1/10 volume of H2O or AcP [freshly prepared in H2O, potas- 
sium lithium salt (Sigma)] at either 100 mM or 1 M concentration 
for a final concentration of 10 or 100 mM AcP, respectively. Reac- 
tions were incubated at 37°C for 90 min. BSA was mixed with 4x 
nuPAGE SDS-polyacrylamide loading buffer (Life Technologies, 
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Naerum, Denmark) and separated by SDS-polyacrylamide gel 
electrophoresis using NuPAGE gels (Life Technologies]. In-gel tryptic 
digestion and peptide isolation were performed essentially as 
described [Shevchenko et al, 1995). For yeast lysate, the reactions 
were stopped by addition of 5-10 volumes — 20°C acetone and incu- 
bated at — 20°C for 1-2 h. Acetone precipitated protein was resus- 
pended in 8 M urea solution (6 M urea, 2 M thio-urea, 10 mM 
Hepes pH8.0) and the concentration determined by Quick Start 
Bradford assay (BioRad). Equal amounts of protein were mixed, 
diluted to 2 M urea by addition of 3 volume 50 mM Hepes, pH8.5, 
proteolyzed by addition of 1:200 (w/w) trypsin protease (Sigma) for 
16-20 h at room temperature. Peptides were purified by C18 
Sep-Pack classic cartridge (Waters, Waltham, MA, USA). Peptides 
were acidified by addition of trifluoroacetic acid (TEA) to 1% and 
applied to a Sep-Pack column that was pre-equilibrated with 5 ml 
acetonitrile and twice with 5 ml 0.1% TEA. The column was 
washed twice with 5 ml 0.1% TEA and once with 5 ml H2O. 
Peptides were eluted with 3 ml 50% acetonitrile in H2O, 100 ^l of 
lOx IP buffer (500 mM MOPS, pH 7.2, 500 mM NaCl] was added, 
and the acetonitrile was removed by vacuum centrifugation to a 
final volume of ~1 ml. Peptide concentration was determined by 
absorbance at 280 nm using a nanodrop (ThermoEisher Scientific] 
spectrophotometer (~6 mg/ml). An aliquot of 200 ^g peptides was 
set aside for subsequent proteome analysis and the remaining 
peptides were used for acetyllysine enrichment (described below). 

Peptide preparation, acetyllysine enrichment, and 
peptide fractionation 

In-solution protein digestion, peptide purification, and acetyllysine 
peptide enrichment were performed essentially as described (Kim 
et al, 2006; Weinert et al, 2013a). Immuno-enriched peptides were 
eluted from anti-acetyllysine antibody resin using acidified H2O 
[0.2% trifluoroacetic acid (TEA)]. Peptide eluates were loaded 
directly onto a strong cation exchange (SCX) microtip column pre- 
pared as described (Rappsilber et al, 2007; Wisniewski et al, 2009) . 
Peptide eluates were briefly evaporated to remove acetonitrile and 
then loaded onto C18 stage-tips as described (Rappsilber et al, 
2007). Proteome measurements were made by fractionation of total 
peptides by the SCX microcolumn method and analyzed by mass 
spectrometry. 

Quantitative mass spectrometric analysis 

We used stable isotope labeling with amino acids in cell culture 
(SILAC; Ong et al, 2002), to measure changes in protein, lysine acet- 
ylation, and phosphorylation abundance. Peptide fractions were 
analyzed by online nanoflow LC-MS/MS using a Proxeon easy nLC 
system (ThermoEisher Scientific] connected to an LTQ Orbitrap 
Velos (ThermoEisher Scientific] or Q-Exactive (ThermoEisher 
Scientific) mass spectrometer. The LTQ Orbitrap Velos instrument 
was operated under Xcalibur 2.1 (ThermoEisher Scientific) with the 
LTQ Orbitrap Tune Plus Developers Kit version 2.6.0.1042 software 
in the data dependent mode to automatically switch between MS 
and MS/MS acquisition as described (Weinert et al, 2011). The 
Q-Exactive was operated using Xcalibur 2.2 (ThermoEisher Scientific] 
in the data dependent mode to automatically switch between MS and 
MS/MS acquisition as described (Michalski et al, 2011; Kelstrup et al. 



2012). All quantitative MS experiments performed in this study are 
summarized in Supplementary Table SI. The mass spectrometry pro- 
teomics data have been deposited to the ProteomeXchange Consor- 
tium (http://proteomecentral.proteomexchange.org] via the PRIDE 
partner repository (Vizcaino et al, 2013] with the dataset identifier 
PXD000507. 

Peptide identification and computational analysis 

Raw data files were processed using MaxQuant software (developer 
version 1.2.7.1) as described (http://www.maxquant.org/; Cox 
et al, 2011). Parent ion (MS) and fragment (MS2] spectra were 
searched against the Saccharomyces Genome Database (SGD) gen- 
ome release r63, January 5, 2010. The search was performed using 
the integrated Andromeda search engine and both forward and 
reversed (decoy) versions of the databases (Cox et al, 2011). Pep- 
tides were additionally filtered for a minimum posterior error proba- 
bility (PEP) score of 0.01, resulting in data sets with estimated false 
discovery rates that were less than the standard 1% used in most 
proteomic studies (Supplementary Table SI). Quantification of TMT 
mass tags may be inaccurate if additional peptides are co-isolated 
with the targeted peptide ion (Ting et al, 2011). In order to mini- 
mize this effect we restricted the quantification of TMT mass tags to 
MS/MS scans in which the targeted parent ion constituted a mini- 
mum of 90% of the total ion current. Mass recalibration was per- 
formed using high confidence identifications based on a initial "first 
search" using a 20 part per million (ppm) mass tolerance for parent 
ion masses and 20 ppm (HCD) or 0.5 Dalton (CID) for fragment 
ions. Spectra were subsequently searched with a mass tolerance of 
5 ppm for parent ions and 20 ppm (HCD) or 0.5 Dalton (CID) for 
fragment ions, with strict trypsin specificity, and allowing up to two 
missed cleavage sites. Cysteine carbamidomethylation was searched 
as a fixed modification, whereas N-acetyl protein and oxidized 
methionine were searched as variable modifications. Where appro- 
priate, acetyllysine was added as a variable modification. 

Acetyl-CoA assay 

Erozen yeast pellets (-100 ODgoo nm] were deproteinized in 400 |.d 
1 N perchloric acid (PGA) containing 13C2 acetyl-CoA as an internal 
standard (ISTD), mixed with 200 |il 0.4 mm acid washed glass- 
beads and vortexed for 2 min at 4°C. Homogenates were centrifuged 
for 10 min, 10 000 g at 4°G. Supernatant (200 |il) was neutralized 
by adding repeated aliquots (10 |il] of 3 M KHC03 under constant 
vortexing, until bubble evolution ceases (90 \A in total). KC104 was 
pelleted by centrifugation for 5 min, 10 000 g at 4°C. Acetyl-CoA 
from the supernatant was measured as described previously, using 
ISTD for calibration (Magnes et al, 2008). Results were normaUzed 
to OD 600 nm measured after extraction. 

AQUA analysis 

Yeast strains expressing GEP-tagged Pgkl and Eas2 at endogenous 
levels from the chromosomal locus (Huh et al, 2003) were 
obtained from Life Technologies. Protein extracts were prepared 
exactly as above when treating whole cell lysate with AcP. After 
treatment with AcP (or mock-treatment) the GEP-tagged proteins 
were enriched using GEP-trap affinity resin (ChromoTek, Martinsried, 
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Germany) and the proteins eluted by boiling in 2x tris glycine 
sample buffer (Life Technologies). Proteins were resolved on pre- 
cast tris glycine gels (Life Technologies) and peptides recovered by 
standard in-gel digestion (Shevchenko et al, 1996). Heavy-labeled 
AQUA peptide standards (AQUA quant pro; ThermoFisher Scien- 
tific) were mixed with Pgkl and Fas2 peptides at the indicated 
concentrations and analyzed by MS. Untreated protein was ana- 
lyzed in four independent experiments and AcP-treated protein in 
two independent experiments, all with comparable results. 

Data analysis 

Gene Ontology (GO) association and enrichment analysis was 
performed using the Database for Annotation, Visualization and 
Integrated Discovery (DAVID) v6.7 (da Huang et al, 2009). Statisti- 
cal tests were performed using R (http://www.r-project.org/ 
index.html). Box plots were generated using Sparklines for Excel 
(http://sparklines-excel.blogspot.com/). The detection limit for nat- 
urally occurring acetylated peptides used to calculate minimum 
increased acetylation in Fig SG was determined by ranking the 
observed "heavy" labeled peptides by peptide intensity and calcu- 
lating the median intensity of the bottom 10%. This intensity value 
was an order of magnitude higher than the least intense "heavy" 
peptide observed in our experiments, thus the estimates of 
increased acetylation based on this detection hmit are conserva- 
tive. 

Supplementary Information for this article is available online: 
http://msb.embopress.org 
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